arXiv:hep-lat/0112001v2 3 Dec 2001 


|hep-lat/0112001 


ANTISYMMETRIC AND OTHER SUBLEADING CORRECTIONS 
TO SCALING IN THE LOCAL POTENTIAL APPROXIMATION 


M. M. TsypinQ 

Department of Physics and Astronomy, Rutgers University, 
Piscataway, NJ 08854, USA 

and 

Institute for Advanced Study, Einstein Drive, Princeton, NJ 08540, USA 


Abstract 

For systems in the nniversality class of the three-dimensional Ising model we com¬ 
pute the critical exponents in the local potential approximation (LPA), that is, in the 
framework of the Wegner-Houghton equation. We are mostly interested in antisym¬ 
metric corrections to scaling, which are relatively poorly studied. We find the exponent 
for the leading antisymmetric correction to scaling uja ~ 1.691 in the LPA. This high 
value implies that such corrections cannot explain asymmetries observed in some Monte 
Carlo simulations. 
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1 Introduction 


This work is devoted to the study of corrections to scaling in the vicinity of the critical 
point: the endpoint of the line of first order phase transition. 

There are many systems that have a critical point belonging to the universality class 
of the three-dimensional (3D) Ising model. In easy-axis magnetic materials that exhibit 
a phase transition between paramagnetic and ferromagnetic phases (these systems most 
directly correspond to the Ising model) the first order phase transition line is in the 
plane (T, h), where T is the temperature and h is external magnetic field. It is located 
at h = 0 and T < Tc, and ends at the critical point (T = Tc, h = 0). 

The system liquid-gas has the line of first order transition, ending at the critical 
point, in the plane (temperature, pressure). This critical point has been very thor¬ 
oughly studied. A comprehensive review, as well as additional references, can be found 
in m. Experimental data, theoretical considerations P, ^ and Monte Carlo simula¬ 
tions of model systems all indicate that it belongs to the 3D Ising model universality 
class. 

The Standard Model of electroweak interactions was shown to have a similar phase 
transition line in the plane (higgs mass, temperature) ii- Monte Carlo simulations 
provide convincing evidence that the endpoint of this line belongs to the 3D Ising 
universality class |0. 

There are strong arguments in favor of a similar first order phase transition line in 
hot nuclear matter, described by quantum chromodynamics (QCD) at high tempera¬ 
ture, in this case in the plane (chemical potential for baryons, temperature) ii- Its 
endpoint is conjectured to belong to the 3D Ising universality class P). 

Monte Carlo simulations play an important role in elucidating the critical properties 
of these systems, especially of more complex ones, such as electroweak matter and nu¬ 
clear matter at high temperature. These simulations normally have to be performed for 
systems of relatively small sizes, usually no larger than 64^ (and often considerably 
smaller), due to their high complexity and the corresponding high computational load. 
This means that the system is not very deep into the scaling region (roughly speaking, 
one cannot reach the correlation length ^ much larger than the system size L), and 
the interpretation of the Monte Carlo data requires proper account for corrections to 
scaling. 

Here we come to the following interesting point. Unlike the Ising model, phase 
transitions in the liquid-gas system, as well as in electroweak and nuclear matter, 
do not have the exact global Z 2 symmetry that corresponds to changing the sign 
of magnetization in the Ising model, i. e. to simultaneous flipping of all spins. For 
example, in the Ising model (without external field) the probability distribution for the 
total magnetization of the system is always perfectly symmetric, despite corrections to 
scaling, while this is not necessarily so in above-mentioned models. The Z 2 symmetry 
is expected to be dynamically restored in the scaling limit [@], but there can be, in 


1 


addition to the usual corrections to scaling, which are even in the order parameter, 
also corrections that are odd. 

Monte Carlo studies of critical points in such systems, which typically aim at ob¬ 
taining the linear mapping of the vicinity of the critical point in the plane of parame¬ 
ters of the model onto the corresponding vicinity of the critical point of the 3D Ising 
model in the (T, h) plane, have to somehow resolve the issue of deviations from the 
Z 2 symmetry 10- Suffice it to say that such mapping requires determination of the 
magnetization-like (M-like) and energy-like (i?-like) directions in the space of observ¬ 
ables, which would be much easier in the presence of Z 2 symmetry: the requirement 
that the probability distribution of the M-like observable should be symmetric greatly 
aids the search for the corresponding direction in the space of observables |Q. 

In practice, the data produced in Monte Carlo simulations demonstrate quite non- 
negligible deviations from Z 2 symmetry, and an interesting question is whether these 
deviations can be attributed to antisymmetric, i. e. M-odd, corrections to scaling 
Such corrections have attracted much less attention in the literature than the usual 
(M-even) corrections; however, several studies in the framework of the e-expansion 
10, as well as renormalization group have been published. 


While the usual M-even corrections to scaling behave as uj = A/z/ k, 0.8, 

L being the characteristic length in the system, the M-odd corrections to scaling are 
governed by their own exponent uja (we prefer this notation to traditional oj^). The 
e-expansion obtained in 
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UJa — 1 + 


11 685 2 107855 + 103680C 3 2 

—e-e H-+ 0(e ). 

6 324 34992 ^ ^ 


( 1 ) 


At e = 1 this series behaves too poorly (1 -|- 1.83 — 2.11 -|- 6.64 ...) to get conclusive 
results; Fade approximants produce the sequence 2.83, 1.85, 2.32 in orders e, e^, e^, 
respectively, leading to the estimate ua ^ 1-5 |H]. The renormalization group compu¬ 
tation [|^] gives uja = 2.4(5). 

These very high values of ooa imply that M-odd corrections to scaling, going as 
should decay much faster than the leading M-even corrections (T““), and quickly 
become negligible when the scaling limit is approached. However, in the Monte Carlo 
study 0 the asymmetry of the probability distribution of the M-like observable at the 
critical point was found to decay much more slowly with the growing lattice size, more 
like L-o-5. 

Before concluding that this observation requires explanation outside the scope of 
corrections to scaling, one has to be conhdent that there are indeed no M-odd correc¬ 
tions to scaling with small exponents (which could conceivably happen, for example, 
if existing computations of uja somehow missed the leading M-odd correction and 
corresponded instead to a subleading correction). To gain such confidence, we have 
performed a computation of M-odd corrections in the framework of the local potential 
approximation (LPA), i. e. the Wegner-Houghton equation [| 
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The WH equation satisfies our needs nicely, as it reproduces robustly all the relevant 
features of the theory. In this approximation the critical exponent r] is zero. 

0, n, 0, in, 


Our study differs from the existing literature on the WH equation 


I9| , p0| in several significant aspects: (1) we are not aware of any previous study of 
M-odd corrections to scaling in the framework of the WH equation, with the exception 
of unpublished work [0; (2) we do not rely on avoiding the singularity when solving 
the WH equation, and thus check whether the large-field domain plays an important 
role in fixing the values of the critical exponents. 

Our main results are as follows. (1) We confirm the absence of M-odd corrections 
to scaling with small exponents. We get the exponent for the leading M-odd correction 
uja ~ 1.691, which is consistent with |0. (2) The large-field domain is not important 
for fixing the critical exponents, at our level of accuracy. 

We conclude that the asymmetry observed in is not explainable by M-odd 
corrections to scaling. 

The paper is organized as follows. In Sect. 2 we discuss the local potential ap¬ 
proximation and the Wegner-Houghton equation. In Sect. 3 we describe our method 
of finding the fixed point and computing the critical exponents. Sect. 4 contains our 
numerical results and conclusions. 


2 The Wegner-Houghton equation 


In this section we remind the reader the structure and the meaning of the WH equation. 
Accurate derivation can be found in original papers 0, 0 and reviews 0, 0 . 


The starting point is the description of the system in the vicinity of the critical point 
by the theory of the 1-component real scalar field 0 with the bare potential V{(j)): 


Z = 


1 




( 2 ) 

( 3 ) 


We will be eventually interested in the three-dimensional case (d = 3). The effective 
potential Uo{ip) in one-loop approximation reads 


Uo{^) = ViV’) + l [ 

Z J\c 


,|<A ( 2 ir)'' 


A , (y 


(4) 


Let us modify the region of integration in the right hand side from all q below the cutoff 
A to g above certain momentum k and below A. This will produce the fc-dependent 
analog of the effective potential, denoted by Uk- 


Uki^) = + 


k<\q\<A {2 tiY 


d‘‘q , (y 


( 5 ) 
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where we have denoted by the area of a d-dimensional unit sphere, divided by (27r)'^. 
For d = 3, Ks = l/(27r^). Taking derivative over k, subtracting the 9 ?-independent 
quantity —^Kdk^~^ In and replacing V{(p) in the right hand side with Uk{p) (“renor¬ 
malization group improvement”), we obtain 

In [z + . (7) 


This equation describes the evolution of the scale-dependent effective potential Uk{p) 
with the change of scale k. For the purposes of the study of the vicinity of the critical 
point, it is convenient to convert it to the form that includes rescaling of p and Uk- 
We introduce the dimensionless parameter t: k = Ae“*, and rescaled quantities Ut and 
Pt- Uk = k^Ut, p = k'^'^pt, where d = 3 is the dimensionality of space, and is the 
dimension of the held p. In terms of t, JJt and pt eq. 0 reads 


dUt 

dt 



+ y^rf-2-2<i<p 



-\- d ■ Ut — dip ■ Pt 


dpt 


( 8 ) 


Using the canonical dimension of p\ dp = {d — 2)12, dropping all the tildes and t- 
subscripts, and denoting derivatives over p by primes, 

+ U") + d ■ U - . (9) 

As in this approximation Z is a constant that depends neither on p nor on t, one can 
conveniently set Z = 1. Restricting to d = 3, we dually get 

f = + (10) 

This is the form of the WH equation that we will use in the following. (One should keep 
in mind that is does not treat properly the additive constant, that is, the (^-independent 
part of U). 


3 Fixed point and critical exponents 

Equation ([T 0 |) describes the evolution of the effective potential U{p) with the change of 
scale. Let us denote its right hand side by F{p). The hxed point is described by effective 
potential that does not evolve with t, that is, U^{p) such that F{p) = const. Usually 
one finds U^i{p) by numerically solving the differential equation F'{p) = 0 and fixing 
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the solution by requiring that it does not run into singularity at large (p [^, |T^ ^ p!5 


Analysis of the evolution over t of small deviations from U^{(p) produces the critical 
exponents. 

Our approach to computation of and critical exponents is as follows. We 

always approximate U{(p) (both U^{ip) and + perturbations) by polynomials up 

to a certain order n in (p: 

n 

j=i 

To determine U^{ip), we search for parameters 02 , 04 ,... a„ of expansion (PT|) that min¬ 
imize the deviation of F{ip) from constant on a reasonably chosen interval. That is, 
we minimize 

( 12 ) 


a = 


2(Pmax —V^max 

At this stage the odd coefficients of expansion (|TT|) are zero, and minimization involves 
n/2 even coefficients of (pT]), plus Fq. The achievable proximity of F{ip) to constant is 
improving rapidly with the increasing order of the approximation (Fig. |l|). 

The next step is the study of deviations of U{ip) from 17* (ip), and their evolution 
over t. As will be clear below, it is convenient to introduce a separate set of parameters, 
bi = b* + 6bi, i = 1.. .n, to parameterize the deviation of U from t/* (which corresponds 
to b*). Then the evolution of parameters follows from 


+ Sh) = F{b- + Sb,) 


F.. 


(13) 


We have to subtract F* = F{b*) in the right hand side, to account for the approximate 
nature of our treatment of the fixed point. Thus 


dU d „ dF „ 


(14) 


Here dU/dbj and dF/dbj are understood to be taken at the fixed point, b*. Defining 
the scalar product 






2 (p 




(15) 


max —^n 


and introducing matrices 




\dbi 


dbjT 


--(I 




(16) 


we obtain 


A ■ ^6b = B ■ 5b. 
dt 


(17) 
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Thus 


(56 = ( 18 ) 

i=l 

where Aj and h) are, respectively, eigenvalues and eigenvectors of A~^B. Aj are exactly 
the critical exponents we are interested in. 

To simplify the computation and improve its numerical stability, we parameterize 
the deviation of U from t/* so that A^j = Sij, namely, 

n 

(7(v>) = f/.(v>) + E'>jA(»’). (19) 

i=i 


where Pj{(p) are Legendre polynomials Pj normalized on [—(^max, 


PAv) = 


' 2j + l 

2 (^ 




max \ Y^max 




( 20 ) 


Then A* are just the eigenvalues of the matrix B. 


4 Results 


The results for critical exponents are collected in Tables 1-3. We use polynomial 
approximations of order 6 ... 16, for three values of (^max: 0.5, 0.44 and 0.4. We ob¬ 
serve that the higher eigenvalues are almost insensitive to the choice of (^max- This is 
fortunate, and provides additional evidence that the LPA captures correctly the im¬ 
portant part of physics. The usual approach relies on avoiding the singularity at large 
tp 1 ^, [^, 1 ^, and one may wonder to what extent the results are determined by 

the asymptotical properties of the WH equation at large 99 , rather than by properties 
at physically relevant range of p (of order of the position of the minimum of f/*). 

Our values of the critical exponents in the even sector. 


yt ^ 1.45041, w ^ 0.5952, cua ^ 2.838, cug ^ 5.18, 
are in agreement with the most accurate of the previous computations [TO 
u ^ 0.6895, a; 0.5952, Aa -2.8384, A 3 ^ -5.1842. 


( 21 ) 


( 22 ) 


The most interesting part of our results is the set of critical exponents for the odd 
sector: 

yh = 2.5, I/shift = 0.5, uja ~ 1.691, a;A 2 ~ 4.00. (23) 

The appearance of the somewhat unusual exponent ?/shift reflects the fact that, strictly 
speaking, in the case of asymmetric U{p) the renormalization group transformation 


6 













should include, in addition to blocking and rescaling of ip, also the shift of (p. As we 
did not take this into account, this additional exponent emerged. 

The value of uja is consistent with but has higher precision. We are not aware 


of any previous computation of u!a 2 - 

To summarize, we have computed the critical exponents, in (p-even as well as in 
(/9-odd sectors, for the systems in the 3D Ising universality class, in the local potential 
approximation. We show that the previously known values of the critical exponents 
in the (/9-even sector are reproduced, within our accuracy, even without relying upon 
avoiding the singularity in the Wegner-Houghton equation at (/9 —>■ -fcx). The absence 
of slowly-decaying (/9-odd corrections to scaling implies that they are not responsible 
for asymmetries observed in Monte Carlo studies [0. 

I am deeply grateful to Prof. H. Neuberger for his kind hospitality, support, and 
many interesting discussions. I would like to thank the Department of Physics and 
Astronomy, Rutgers University and the School of Natural Sciences, IAS, for their hos¬ 
pitality and support. This work was supported in part by DOE grants DE-FG02- 
96ER40949 and DE-FG02-90ER40542. 
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n 

a 


16 

1.48x10“® 

14 

2.36x10“'^ 

12 

8.48x10“’^ 

10 

2.31x10“® 

8 

4.13x10“® 

6 

2.78x10“"^ 

Vh 

A 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

yt 

S 

1.450413 

1.450412 

1.450402 

1.45044 

1.44878 

1.44099 

2/shift 

A 

0.5 

0.500003 

0.500035 

0.500223 

0.505546 

0.491122 

—UJ 

S 

-0.595235 

-0.595256 

-0.595119 

-0.595044 

-0.589349 

-0.530546 

-^A 

A 

-1.69132 

-1.6914 

-1.69099 

-1.70327 

-1.5294 

-0.920653 

—102 

S 

-2.83838 

-2.83869 

-2.83794 

-2.89125 

-2.50668 

-2.00335 

—^A2 

A 

-3.99828 

-4.00259 

-4.11481 

-3.46038 

-3.38734 


—(03 

S 

-5.18418 

-5.20558 

-5.37167 

-4.48014 

-5.0081 



A 

-6.44866 

-6.64226 

-5.61029 

-6.89618 




S 

-7.76251 

-7.95578 

-6.91652 

-9.08612 




A 

-9.35137 

-8.46459 

-11.5582 





S 

-10.8738 

-10.2681 

-14.3183 





A 

-12.3235 

-17.3682 






S 

-14.6305 

-20.7111 






A 

-24.3498 







S 

-28.2887 







Table 1: Critical exponents for (prasa = 0.5. The qnality of approximation is character¬ 
ized by a, according to (0)- Second column indicates the symmetry of the eigenvector 
{S = even, A = odd). 
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n 

a 


16 

3.41x10“® 

14 

3.96x10“® 

12 

3.14x10“'^ 

10 

3.43x10“'^ 

8 

1 . 41 x 10“5 

6 

1.28x10“"^ 

Vh 

A 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

yt 

S 

1.450410 

1.45035 

1.45056 

1.4505 

1.44882 

1.45284 

I/shift 

A 

0.499991 

0.499918 

0.500252 

0.500054 

0.498699 

0.523208 

—OJ 

S 

-0.595165 

-0.595677 

-0.594694 

-0.593543 

-0.604152 

-0.546673 

—UJa 

A 

-1.69158 

-1.69361 

-1.68492 

-1.7014 

-1.75849 

-1.92471 

—U2 

S 

-2.83969 

-2.84017 

-2.83152 

-2.90894 

-3.08336 

-3.59606 

—^A2 

A 

-3.99407 

-4.03145 

-4.2279 

-4.56781 

-5.65287 


—UJs 

S 

-5.19053 

-5.32954 

-5.73614 

-6.33275 

-8.01733 



A 

-6.77786 

-7.47611 

-8.36842 

-10.7213 




S 

-8.44641 

-9.46912 

-10.6988 

-13.8132 




A 

-11.7329 

-13.3336 

-17.2785 





S 

-14.2826 

-16.2727 

-21.1253 





A 

-19.5243 

-25.361 






S 

-23.0967 

-29.9828 






A 

-34.9966 







S 

-40.4086 







Table 2: Critical exponents for v^max = 0.44. 
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n 

a 


16 

1.65x10-® 

14 

6.53x10“® 

12 

1.15x10“'^ 

10 

3.14x10“'^ 

8 

5.78x10“® 

6 

6.64x10“® 

Vh 

A 

2.5 

2.5 

2.5 

2.5 

2.5 

2.5 

Vt 

S 

1.45031 

1.45014 

1.45136 

1.45019 

1.44437 

1.47603 

I/shift 

A 

0.499857 

0.499766 

0.501045 

0.499452 

0.493915 

0.537914 

—ui 

S 

-0.595091 

-0.597076 

-0.593018 

-0.590154 

-0.619136 

-0.6057 

—UJa 

A 

-1.69522 

-1.6978 

-1.66645 

-1.73619 

-2.07471 

-2.76637 

—U2 

S 

-2.84795 

-2.83693 

-2.82088 

-3.05967 

-3.72532 

-4.91745 

—^A2 

A 

-3.98292 

-4.09967 

-4.59586 

-5.64517 

-7.68927 


—Co’s 

S 

-5.20884 

-5.56624 

-6.46 

-7.98584 

-10.7612 



A 

-7.30803 

-8.68827 

-10.7014 

-14.2306 




S 

-9.3943 

-11.2751 

-13.8281 

-18.2123 




A 

-14.2416 

-17.3769 

-22.65 





S 

-17.6035 

-21.3259 

-27.5739 





A 

-25.6903 

-32.9993 






S 

-30.4805 

-38.8901 






A 

-45.2637 







S 

-52.1331 







Table 3: Critical exponents for v^max = 0.4. 
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U* (tp) 





Figure 1: Approximate fixed point effective potential U^{ip) and the corresponding 
deviation of the right hand side of eq. m from constant, for n = 6 and n = 16 
(‘^max = 0.5 in both cases; the two graphs for U^{(p) are hardly distinguishable). 
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